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Formation and fragmentation of networks is typically studied using percolation theory, but most previous 
research has been restricted to studying a phase transition in cluster size, examining the emergence of a giant 
component. This approach does not study the effects of evolving network structure on dynamics that occur at 
the nodes, such as the synchronization of oscillators and the spread of information, epidemics, and neuronal 
excitations. We introduce and analyze new link-formation rules, called Social Climber (SC) attachment, that 
may be combined with arbitrary percolation models to produce a previously unstudied phase transition using the 
largest eigenvalue of the network adjacency matrix as the order parameter. This eigenvalue is significant in the 
analyses of many network-coupled dynamical systems in which it measures the quality of global coupling and is 
hence a natural measure of connectivity. We highlight the important self-organized properties of SC attachment 
and discuss implications for controlling dynamics on networks. 

PACS numbers: 64.60.ah, 89.75.-k, 87.23.Ge 



I. INTRODUCTION 

Dynamics on networks has become a research area of broad 
importance, with considerable effort focused on understand- 
ing how dynamics are affected by network structure lHl-ll7l. 
Of particular interest are dynamics that depend on global mea- 
sures of network connectivity, and in particular on the largest 
eigenvalue A of the network adjacency matrix A (Aij ^ if a 
link exists from node i to node j). We will refer to this broad 
class, which includes models for synchronization [1], genetic 
expression J2l, neural excitation H, and epidemic spreading 
10. IH, as connectivity-governed dynamics. We note, however, 
that while analyses of such systems JU-Hl] typically assume 
that the network structure is static and connected (i.e., lack- 
ing isolated nodes/clusters), many applications exist for which 
the network structure is non-static and/or fragmented, such as 
epidemic spreading with immunization |5l|6[], communication 
and transit systems operating under failure or attack [0], and 
information processing in the brain flU. 

The systems that we categorize as having connectivity- 
governed dynamics share a common property of being eas- 
ily manipulated through changes in A. For example, one can 
prevent viral spreading in technological and social networks 
by decreasing A through immunization 15|] or promote the dis- 
semination of information in communication and sensor net- 
works by increasing A Hfl. In recent years there has been 
much interest in studying the effects of topological modifi- 
cation on A and developing efficient strategies for tuning A 
through the addition and/or subtraction of links and/or nodes 
H-Qj]]. However, such perturbation techniques do not address 
networks undergoing formation or fragmentation processes, a 
prob lem traditionally studied with network percolation theory 
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From network resilience against targeted attacks and fail- 
ures lfl3l [Till to dynamics on networks which often 
map spreading dynamics to percolation processes, there are 
widespread applications for percolation theory in the field of 
complex networks. However, applying such techniques to dy- 
namics on evolving networks is hindered in that it can require 
two levels of analysis, theory for the change in network topol- 
ogy and theory for the dynamics. We hypothesize that a cen- 
tral element frustrating the development of this approach is 
that this field of research has largely focused on studying a 
phase transition in the size of the largest connected cluster by 
analyzing cluster aggregation and the emergence of a giant 
component, a cluster whose size is of the same order as the 
entire system lfl2ll . However, the application of subsequent 
percolation theory requires information about cluster topol- 
ogy in addition to cluster size, which highlights the need for 
percolation theory focusing on other cluster properties such 
as the spectra of clusters (i.e., the eigenvalues of adjacency 
matrices corresponding to clusters), modularity, assortativity, 
transitivity, etc. lfl5ll . 

Here we study the link-percolation phase transition using A 
as a novel order parameter, shedding light on a new phase tran- 
sition in connectivity, corresponding to a poorly-connected 
network becoming well-connected (or vice versa) in terms of 
the topology's effect on dynamics. In order to produce such 
a transition, we introduce a link-formation rule called Social 
Climber (SC) attachment for which we derive the asymptotic 
scaling behavior of A for large network size N. We show that 
networks forming under SC attachment exhibit maximal scal- 
ing A ~ 0(N 1 ^ 2 ), indicating that our model may be of broad 
interest for the design of networks with large A, a property that 
is often beneficial |9j] and can lead to, for example, excellent 
robustness against attack and failure lfl3l [Till and very good 
spreading characteristics lUl 0]. It follows that SC attach- 
ment is a promising approach for the design of self-organized 
communication and sensor networks IT6ll with topologies de- 
signed for the rapid dissemination of information. We demon- 
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strate this application by showing that networks formed un- 
der SC attachment exhibit enhanced spreading properties with 
respect to the susceptible-infected-susceptible (SIS) model 
m5[], a contagion model with many applications including the 
dissemination of information, sometimes referred to as "gos- 
sip based" communication or epidemic routing lfl7ll . We note 
that the development this potential application may be facili- 
tated by the fact that SC attachment may be combined with ar- 
bitrary percolation processes, such as Erdos-Renyi (ER) per- 
colation lfl8ll and Achlioptas processes lfl9l - l2lll . to indepen- 
dently control cluster aggregation (determined by the percola- 
tion process) and connectivity within clusters (determined by 
SC attachment). 



II. SOCIAL CLIMBER ATTACHMENT 



SC attachment DSC attachment 




(a) (b) 

FIG. 1 . (color online) Link x is proposed by an arbitrary percolation 
model (dashed lines) to connect nodes a and b, merging clusters C a 
and C'b- Because C a 7^ Cb the proposed link x is discarded and 
instead one of the following new links is formed, (a) Node a is linked 
to the largest-degree node in C'b with link y to model SC attachment. 

(b) The largest-degree nodes from C' a and Cb are linked together with 
link z to model DSC attachment. These processes may be visualized 
using free PercoVIS software I22I1 . 



A link-percolation process begins with N isolated nodes, 
indexed n = 1, 2, . . . , N. In discrete steps t = 1, 2, . . . , a 
new undirected link between two nodes is selected according 
to a rule or set of rules, and is then formed. Thus, after i 
steps there will be i links in the network, resulting in clusters 
of connected nodes, each of whose size (number of nodes in 
the cluster) may range from one (an isolated node) to N (a 
cluster that spans the entire network). Depending on the rules 
used to select links, the evolution of cluster sizes, and in par- 
ticular the size of the largest cluster, may vary significantly. 
Social Climber attachment introduces a new link reselection 
step between link selection and link formation, which we mo- 
tivate by analogy to a corresponding social process: collo- 
quially, a "social climber" is someone who actively attempts 
to make powerful friends in order to become more powerful 
himself. When introduced to a new person, a social climber 
learns about the relative popularity of the people in that per- 
son's clique and eventually befriends whoever is of maximal 
importance. With this in mind, SC attachment is a link res- 
election step during percolation where the proposed link be- 
tween two nodes is altered by allowing one of those nodes to 
act like a social climber, choosing to link to the node of max- 
imal importance in the other node's cluster. Therefore, given 
a link-percolation process, we summarize SC attachment as 
follows, (i) Let x be a proposed undirected link connecting 
nodes a and b, generated by an arbitrary percolation model, 
(ii) Let clusters C a and Cb be the clusters to which nodes a 
and b belong, respectively. Then, if C a 7^ C& the proposed 
link x is discarded and instead a link y is made between node 
a and the largest-degree node in Cb, as shown in Fig.QJ. (iii) 
If nodes a and b belong to the same cluster, C a = Cb, then 
the proposed link x is made without modification. Note that 
SC attachment does not affect which clusters combine, but 
does affect the topology of the resulting joined cluster. The 
SC model chooses a connection to the node of largest degree 
in a cluster using nodal degree as a proxy for the dynamical 
importance measure, DI = u n v n iflOfl ^ where u and v denote 
the right and left eigenvectors of A corresponding to A. For 
the undirected networks considered here, symmetry of A im- 
plies u = v, and thus the node with largest eigenvector entry 
u n will have maximal DI in its cluster iflOfl. Provided that the 



node with largest eigenvector entry u n also has largest degree 
k n , we allow the SC model to select nodes based on degree 
for simplicity and ease of computation. One may equivalently 
view SC attachment as forming a link to the node with largest 
degree by using degree centrality as a proxy for eigenvector 
centrality. 

In addition to SC attachment, we introduce Double Social 
Climber (DSC) attachment, in which a proposed link between 
nodes a and b is either replaced by a link between the nodes 
with maximal degree in each cluster when C a 7^ Cb, as shown 
in Fig. [Tp, or is formed between a and b without modification 
when C a — Cb- We note that DSC attachment corresponds to 
maximizing connectivity of the resulting cluster, as measured 
by A, whenever the node of maximal degree is also the node 
of maximal eigenvector entry in each cluster SC and DSC 
attachment may be visualized using free PercoVIS software 



in. ANALYSIS 

Although we will later generalize our methods to other per- 
colation rules, we first analyze SC and DSC attachment for the 
well-known ER percolation process 01811 . The rule for select- 
ing a link in ER percolation is simple: two nodes are chosen 
uniformly at random and a link is formed if there is not al- 
ready a link between them. Traditional analysis has focused 
on the relationship between the number of links added t and 
the size of the largest cluster G(t), called the giant component 
(GC) when G{t) ~ 0{N). It is convention to rescale both i 
and G by N [i.e. t = i/N and G{t) = Gti)/N], where one 
obtains in the asymptotic limit N — s- 00 il8fl . 

G(t) = { 1 - e- 2tG W ] \ > ol ■ (1) 

Here, for variable control parameter t, the network undergoes 
a second order phase transition in cluster size at the perco- 
lation threshold tf B = 0.5, as observed through the order 
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parameter G(t). Because SC attachment affects the topology 
of clusters and not their sizes, Eq. ([1} remains valid for ER 
percolation combined with SC attachment. 

We begin our analysis by studying the emergence of large- 
degree nodes. For a given time t, consider a large cluster C 
containing s 3> 1 nodes, and let fc max denote the maximal 
nodal degree in C. By large cluster, we mean that s is with 
high probability larger than the size of another randomly cho- 
sen cluster in the network, and eventually we will consider 
only the case in which cluster C is the largest cluster in the 
entire network. We will compute the expected change in /c max 
for the addition of a single link. When a link is proposed be- 
tween nodes a and b by ER percolation, k mdx will increase 
by one if: (i) a ^ C and b € C (depicted in Fig. QJ where 
Cb = C), or (ii) a,b g C and the degree of a or b is fc max . 
Since ER percolation chooses nodes uniformly at random, the 
probability that a^Cisl — s/TV, and the probability that 
b € C is s/N. Since a and b are chosen independently, the 
probability of case (i) is (s /TV) [1 — s/N]. The probability that 
a randomly chosen node in a cluster of size s has degree k mdx 
is r/s, where r is the number of nodes in that cluster with de- 
gree fc max . Thus, the probability of case (ii) is, to leading order 
as TV -> oo, (s/N)[l-s/N](r/s) + (s/N) 2 [2r/s], where the 
first term corresponds to a 6 C, b ^ C and the second term 
corresponds to a, b G C. We note that other corrections may 
be included to address the chance in (i) that the maximal de- 
gree of C a is larger than the maximal degree of Cb, k max + 1, 
but such corrections decay rapidly as the size difference be- 
tween Cb and C a increases. Since s 3> 1 while r ~ 0(1), 
case (i) is the dominating process for sufficiently large s, so 
the expected rate of change of fc max averaged over all possible 
links is 
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Here, £ = 1 for SC attachment and £ = 2 for DSC attachment, 
since for DSC case (i) applies to both a ^ C, b € C and 
a G C, b £ C. Using d/di = Nd/dt, integration of Eq. © 
predicts that the largest degree of a node within the GC at time 
t is 



CN / [G(r) - G 2 (r)] dr. 



(3) 



Note that this scaling with TV is the largest achievable scaling 
of a degree. For comparison, in networks with power-law de- 
gree distribution, P(k) oc k~~< , the expected maximal degree 
scales as ©(TV 1 /^ -1 )), approaching 0(N) as 7 -> 2+. 

In o rder t o understand the implications of Eq. (|3), we use 
A ss Vfc max , an asymptotically (N — > 00) accurate approxi- 
mation derived in E4ll and discussed further in |)25]] . While the 
model used to generate this estimate is not equivalent to SC 
attachment, we find it remains accurate here. Using this esti- 
mate in conjunction with Eq. (0, and noting that for t > tf R , 
the largest eigenvalue of the GC will be larger than the largest 
eigenvalues of smaller clusters, we obtain the following ex- 
pression for the expected largest eigenvalue for ER percola- 




FIG. 2. (color online) Xer is shown for ER percolation with DSC 
attachment. Our prediction, Eq. (|4j with G(t) given by Eq. ([]} 
(solid line), agrees well with observed values for simulations with 
N = 10 5 (X symbols) and N = 10 6 (circles). Observed values 
for Xae (crosses) also agree well with theory, Eq. ((6) (dashed line), 
for AE percolation with DSC attachment for TV = 10 6 . (inset) For 
TV = 10°, the great extent to which SC and DSC attachment increase 
connectivity is shown by comparing to A for classical ER percolation. 



tion with SC attachment, 



E[\ ER (t)} = \ (N / [G(T)-G^r)}dr, 



(4) 



implying that the network undergoes a continuous phase tran- 
sition in connectivity at precisely the same value t = tf R at 
which a phase transition in cluster size occurs. In the super- 
critical regime, A achieves maximal scaling, A ~ 0{N 1 ' 2 ). 
For comparison an all-to-all network has similar scaling, 
A = y/N - 1, but uses 0{N 2 ) links compared to 0(N) 
used by SC attachment. Asymptotic scaling constants of 
Eq. © for large t may be solved by integrating with re- 
spect to G, rather than t, and using a dilogarithm to obtain 

\ ER (t)/VN 



7r 2 /12 w 0.42 for SC attachment and 
^2-71-2/6 « 0.6 for DSC attachment 
We confirm the accuracy of Eq. by direct simulation of our 
model, shown in Fig. [2] which demonstrates excellent agree- 
ment between Eq. (HJi (solid line) and observed values of Xer 
for ER percolation with DSC attachment with TV = 10 5 (X 
symbols) and TV = 10 6 (circles). For integration in Eq. (3) we 
use the asymptotic theoretical value G(t) given by Eq. (Q~|). In 
the inset of Fig. [2] we compare observed values of Xer for SC 
and DSC attachment to traditional ER percolation, where our 
models' main effects are highlighted: under SC or DSC at- 
tachment Xer attains significantly larger values and also un- 
dergoes a sharp increase at tf R . 

Because the expected connectivity, measured by Xer, is 
a function of G(-), developing scaling arguments for Xer is 
straightforward since scaling for G(t) is known: G ~ log(TV) 
for t < tf R and G - TV for t > if R fH. Therefore, when 
SC is used in conjunction with ER percolation, in the limit 
TV — > 00 we have 



\ ER {t) 

\E R (t) 



TV 1 / 2 
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(logTV) x l 2 scaling 
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FIG. 3. (color online) Scaling of \Enit) predicted in Eq. ((5} is 
demonstrated using the ratios of \ E B.{t) for TV = 10 6 and TV = 10 5 , 
denoted <f>Eii(t). Agreement between the prediction of Eq. ((5} (lines) 
and measurement from simulation (crosses) is good. Measurements 
from AE percolation, 4>AE(t), are also shown (circles). 



To validate Eq. (O we estimate the change in X E r when 
system size is increased by defining <fi E R.(t) as me rat i° °f 
X ER (t) for TV = 10 6 to X E R(t) for TV = 10 5 . As shown in 
Fig. [3j we predict that <f> E R(t) « y^6/5 for t < if R (dashed 
line) and 4> ER (t) for t > tf R (solid line), both of 

which agree well with <t> E R,{t) calculated from a single simu- 
lation of each system size (X symbols). 

The methods used to derive Eqs. (|3]|5j, which involved cal- 
culating the probability that an isolated cluster attaches to the 
GC, may be easily adapted to other percolation models. For 
example, consider Achlioptas processes 1IT9I - I21I1 for which 
the merging of clusters depends on cluster size (up to some 
bound). This class of percolation models has recently re- 
ceived much attention, focusing on analysis of a rapid phase 
transition in cluster size referred to as "explosive percolation" 
lEol EH l23ll . Repeating the reasoning process in deriving 
Eqs. d3j5j for Adjacent-Edge (AE) percolation 1I20I1 , we pre- 
dict the largest eigenvalue to be 



E[X AE (t)} =^J* [G(t) + G*(t) - 2G3(t)] dr, (6) 

where again £ = 1 for SC and C = 2 for DSC. Note that 
despite the near-discontinuous phase transition in G, maximal 
scaling A ~ ©(TV 1 / 2 ) is still achieved. In Fig. [2] we show 
good agreement between observed values for Xa E (crosses) 
and Eq. © for DSC (dashed line), where observed values for 
G(t) were used in Eq. © as an analytic expression has yet 
to be developed. Note that X AE (t) < X E n(t), which we at- 
tribute to the integrands of Eqs. (|4]i and which are max- 
imized at G = 1/2 and G = (1 + V<)/6 respectively and 
are zero at G = and G = 1. Since AE percolation pro- 
duces rapid growth in G, the integrand of Eq. (|6]l is not large 
over a majority of the integration interval, so essentially, the 
explosive growth in G minimizes the regime during which SC 
attachment has a large effect on A. In Fig. [3] we also show 
4> AE {t), the ratio of X AE (t) for TV = 10 6 and TV = 10 5 , 
where we observe similar scaling in the subcritical and super- 
critical regimes as observed and predicted for X E r. 



IV. EXPERIMENTATION 

We now demonstrate the effect of SC attachment on dynam- 
ics. Because SC attachment produces networks with maximal 
scaling of A, we focus on an application in which large A is 
beneficial: the dissemination of information in communica- 
tion and wireless sensor networks [[l6||, which is often mod- 
eled as an epidemic |[T7ll . We note, however, that large A is 
not always advantageous. For example, large A in ecolo gica l 
networks can promote instability and species extinction 112711 . 
(See U for a discussion of applications in which it is benefi- 
cial to have either small or large A.) Here we study SIS con- 
tagion J^Ht], which has been used to study spreading process 
from viral propagation in social and technological networks 
to the dissemination of information such as rumors and data 
d. 

To briefly review, the SIS model is a continuous time pro- 
cess in which each node may be susceptible to infection or 
infected. Each infected node may infect each of its suscep- 
tible network neighbors at rate a, and each infected node 
may also spontaneously heal and return to being susceptible 
at rate /3. The network state in which no nodes are infected 
and all nodes are susceptible is a fixed point of the collective 
dynamics, but this fixed point may not be stable to perturba- 
tion (i.e., a small fraction of nodes being infected by some ex- 
ternal agent). For many topologies of connected networks in 
which a fraction of nodes are initially infected, the expected 
steady-state fraction of infected nodes / may either be zero 
(no infections, stable fixed point) or nonzero (endemic infec- 
tion, unstable fixed point), depending on whether a/ (3 sur- 
passes the epidemic threshold A -1 14j,|5[]. Note that endemic 
infection can be prevented by decreasing A through immu- 
nization until A -1 > a//3. Interestingly, for very large infec- 
tion rates, reducing A to prevent endemic infection can require 
the complete fragmentation of the network. For example, if 
a/(3 > 0.5 the prevention of endemic infection requires A < 2 
which guarantees fragmentation of the network IU4II . This sce- 
nario has been observed experimentally for virus propagation 
on mobile phone devices 12811 , where slowly spreading Blue- 
tooth viruses may be inhibited by immunization (i.e., antiviral 
software) but rapidly spreading messaging viruses are inhib- 
ited only by a fragmented network. 

We simulated SIS dynamics for moderate a/ {3 on two net- 
works forming under ER percolation, one with SC attachment 
and the other without, predicting that SC attachment will have 
a significant impact on the steady state fraction of infected 
nodes, /. We simulated dynamics with (a,/3) = (0.075,1) 
on TV = 10 5 nodes at many points t in the percolation process, 
initially infecting 1 % of nodes and then allowing the system 
to reach a steady state fraction of infected nodes f(t), before 
allowing percolation to continue to another value of t, where 
the dynamics were re-initialized, simulated, and so on. The 
resulting curves f(t), are shown in Fig. |4] where the shaded 
region highlights that networks forming with SC attachment 
(open squares) have significantly enhanced spreading charac- 
teristics compared with networks forming without SC attach- 
ment (filled squares). To contrast this result, we also plot f(t) 
for a large infection rate (a,/3) = (0.5, 1), where f(t) with 
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FIG. 4. (color online) SIS epidemics |Q] were simulated on a net- 
work forming by ER percolation with SC attachment (open symbols) 
or without (filled symbols), where tN is the number of links added. 
f(t) is the steady-state fraction of infected nodes when 1% of nodes 
are initially infected. The shaded region highlights the significant 
impact of SC attachment for moderate infection rate (squares). For 
high infection rate, SC attachment has no effect (circles). 



SC attachment (open circles) is indistinguishable from that 
without SC attachment (filled circles). This is not surprising 
as one would expect any initial infection to saturate the cluster 
in which it begins, in which case f(t) would depend primarily 
on cluster size, not topology. We thus find two regimes of SIS 
dynamics on fragmented networks: when a/ (3 is sufficiently 
large, / depends primarily on network fragmentation (i.e. the 
size of the GC), but for moderate and small a/ j3, f depends 
strongly on the connectivity of clusters (i.e., their respective A 
values). 



V. DISCUSSION 

Motivated by the need for the development of analysis for 
connectivity-governed dynamics OrlH] on evolving networks, 
we have developed a percolation theory focusing on the con- 
nectivity of clusters, rather than their size. In this pursuit, we 
have introduced a model, Social Climber attachment, that pro- 
duces networks with strong connectivity and maximal scal- 
ing of A, and validated our claims using two link-percolation 
models. While strong connectivity in networks is achievable 
via other percolation models (e.g., networks with heavy-tailed 
degree distributions generated by the Chung-Lu model ll29ll ). 
such methods typically require that the nodal degrees and net- 
work connectivity be defined a priori. In contrast, two key 
properties distinguish SC attachment, (i) First, SC attachment 
produces networks with large A via self-organization. Be- 
cause A governs many dynamical processes ilj-|5|], SC attach- 
ment provides a foundation for designing networks that self- 



organize with properties linked to large A such as robustness 
lfl3i[T4ll and the efficient spread of information H. Our model 
is therefore promising as a starting point for the development 
of self-organized communication networks such as wireless 
sensor networks 111 611 . where data broadcasting may be mod- 
eled by SIS transmission lfl7tl . Development and analysis of 
this application may be facilitated by the fact that SC attach- 
ment does not affect cluster sizes, only their internal topology, 
(ii) Second, a novel phase transition in connectivity occurs 
for networks forming under SC attachment, which may have 
broad applications. For example, A(i) changes most rapidly 
near the percolation threshold, so creating networks near criti- 
cality may offer an effective approach for designing networks 
on which dynamics can be efficiently controlled by adding 
or removing a minimal number of links. This approach may 
therefore aid in the design of critical infrastructure (e.g., the 
power grid, communication networks, and airline networks) 
that can be easily switched between topologies designed for 
high-flow and low-flow conditions. 

We conclude by suggesting several possible extensions to 
this work that may be of interest to readers. First, SC attach- 
ment uses complete information about the structures of the 
clusters that it connects, yet in some applications this infor- 
mation may be difficult or impossible to obtain. The effects 
of incomplete information or noise on SC attachment are as 
yet unexplored. For example, incorporating a probabilistic 
(rather than deterministic) link reselection step may be of in- 
terest as the resulting process would have some similarity with 
the preferential attachment network growth model [|30ll . Sec- 
ond, one may wish to adapt our model to study various real- 
world networks. While networks formed under SC attachment 
feature large A, to incorporate this model for the design of 
communication networks one would likely need to consider 
many other design criteria such as betweenness centrality and 
software protocols lfl6ll . Finally, we named our model Social 
Climber attachment to reflect the selfish behavior of individu- 
als in social situations, yet the generation of a network topol- 
ogy similar to that observed in social networks using the SC 
model would require additional link-formation rules, such as 
those producing modularity and transitivity lfl5l . 
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